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Abstract 

°. 

We have modified the Herwig event generator to incorporate diffractive 
' interactions. All standard Herwig hard subprocesses are available. 

j3 \ 1 Introduction 

> 

"~ 1 ' rim 

^ The modifications to Herwig [flj that are necessary to incorporate diffractive 

collisions are relatively simple once it is noticed that pomeron exchange events in 
hadron - hadron collisions look very much like resolved photoproduction events 
in lepton - hadron collisions. In resolved photoproduction in electron - proton 
collisions, for example, the process is modeled by the incoming electron radiating 
a quasi-real photon according to a flux formula. The photon is then treated as 
a hadronic object with a structure function, which undergoes a collision with the 
beam proton. Similarly, single diffractive interactions in proton - proton collisions 
may be modeled by assuming that one of the beam protons emits a pomeron, again 
according to some flux formula, which subsequently undergoes an interaction with 
the other beam proton (see figure 1). Herwig will automatically choose to radiate 
a photon from a beam lepton if a hard subprocess is selected which requires a 
hadronic structure for the beam lepton. Examples would be choosing Herwig 
subprocess 1500 (QCD 2 — > 2 scattering) in an electron - proton collision, or 
subprocess 9000 (deep inelastic scattering) in an e + e~ collision. 

All that is necessary therefore, to simulate a single diffractive interaction in a pp 
collision is to replace the photon flux with a suitable pomeron flux factor, and the 
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Figure 1: The POMWIG philosophy. Photoproduction in ep (or ee) collisions is 
replaced with pomeron or reggeon exchange in pp (or ep) collisions. 



photon structure function with a pomeron structure function, and run HERWIG 
in ep mode rather than pp mode. The electron should then be identified with 
the proton which remains intact after the diffractive scattering. This process may 
be generalised to include sub-leading exchanges, and to perform double pomeron 
collisions and diffractive DIS as will be described below. 

Our philosophy for Pomwig has been to make as few changes as possible to the 
Herwig code. In particular, the Herwig common blocks are left unchanged, 
and only two Herwig subroutines have been modified. This approach has some 
disadvantages, in that we have left the event record unchanged. An intact proton 
will still appear as a lepton, therefore, and a pomeron will appear as a photon. We 
felt that such cosmetic inconveniences are outweighed by the overall simplicity and 
ease of installation and maintenance which are features of Pomwig . 



2 Installing the Code 

The code can be obtained on request from the authors, or from ||. The routines 
supplied will function with all currently available versions of Herwig from 5.9 
onwards. There are two HERWIG routines which must be replaced: HWSFUN and 
HWEGAM. HWSFUN has been modified to call the HI diffractive structure function 
routine H1QCD ||, or the user defined structure function routine POMSTR in the case 
of pomeron exchange, or to use a neutral pion structure function in the case of 
reggeon exchange. The pion structure is taken from the CERN pdf lib package 
via a call from HWSFUN. The default reggeon structure is that of Owens [|J, as 
used in the HI diffractive DIS analysis 0. HWEGAM substitutes the HI pomeron or 
reggeon flux for the photon flux. The flux itself is calculated by the new routine 
FLUX. 

If beam particle 1 is the electron (playing the role of the diffracted proton), then 
setting HERWIG parameter MODPDF(l) = -1 and NSTRU = 6 will cause the HI 
pomeron structure function to be used for beam particle 1. Similarly, setting 
NSTRU = 7 will select reggeon exchange, and NSTRU = 8 will select the user defined 
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pomeron structure from POMSTR. The parameters in the HI pomeron structure 
function routine are set via an initialisation call to H1QCD. 

To initialise the HI pomeron structure function, the call should be made from the 
HERWIG routine HWABEG as follows; 

INTEGER if it 

DOUBLE PRECISION XPQ(-6:6) ,X,Q2 

02=75 

if it=5 

X=0.1 

call QCD_1994(X,Q2,XPQ,if it) 

Details of the parameters can be found in the H1QCD routine. The above parameters 
select the LO fit 2. For this initialisation call to H1QCD, the values of Q2 and X are 
irrelevant. 

The HERWIG parameters 02WWMN and 02WWMX set the minimum and maximum 
pomeron virtualities, t m in and t max respectively. YWWMIN and YWWMAX set the min- 
imum and maximum incoming proton energy fractions carried by the pomeron, 
and x^ ax respectively. 

The default parameters for the pomeron and reggeon fluxes are those found by the 
HI Collaboration in ||, for the case in which no interference is assumed between 
the pomeron and reggeon contributions to F^ 3 \ as shown in table l[]. The fluxes 
are parameterised as 

fp/ P (xjp) = N (1) 

J tm ax Xjp 

rt m in e Pm(t) 

fn/p(xp) = C M 2a R (t)-l ( 2 ) 

J tmax Xjp 

where ajp(t) = op(0) + a' P t and or(£) = ojr(0) + a M t. In the case of the user 
defined structure functions, N = 1 by default, and the normalisation may be 
included within POMSTR or FLUX at the users discretion. The normalisation of the 
flux is arbitrary in the case of the HI pomeron structure function. The H10CD 
routine is implemented such that the generated cross section will always match 
F^ 3 ^ as measured by HI at xjp = 0.003, irrespective of the parameters chosen for 
the flux. 



■"■Note that the reggeon parameters are not well constrained by current measurements. The 
defaults are those used to generate the curves in figure 2. 
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Quantity 


Value 


OLp 


1.20 


Ot-IR 


0.57 


ajp 


0.26 


«R 


0.90 


Bjp 


4.6 


Bm 


2.0 


Cm 


48 



Table 1: The default parameters in Pom wig 



Process 


BEAM1 


BEAM2 


MODPDF(l) 


MODPDF(2) 


NSTRU 


PP -> PX 


P 


£+ 


pdflib 


-1 


6 (7, 8) 


PP -» PPX 


E+ 


P+ 


-1 


-1 


6 (7, 8) 


eP -» ePA" 


E+ 


E+ 


N/A 


-1 


6 (7, 8) 



Table 2: The diffractive processes available in Pomwig 



3 Summary of available processes 

The processes currently available in Pomwig are summarised in table 2. The 
particular hard subprocess required is set as usual via Herwig parameter iproc. 
Note that at present for double pomeron exchange, it is only possible to simulate 
pomeron - pomeron or reggeon - reggeon collisions, and not pomeron - reggeon. 
In diffractive DIS (iproc= 19000), the pomeron is always emitted from the second 
beam particle. Note also that process eP — > ePX in table 2 refers only to diffractive 
DIS. Pomwig is unable at present to generate this process with a resolved photon 
component. 

4 Example processes 

In order to make the modifications clear, we shall consider the simulation of several 
example processes in both ep and pp collisions. First, we consider the diffractive 
structure function as measured by the HI Collaboration at HERA ||. The 

complete list of modifications to the HERWIG main program in order to generate 
pomeron exchange events using the HI pomeron structure are as follows : 
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PROGRAM HWIGPR 
C— BEAM PARTICLES 
PART1= , E+ ' 

PART2='E- ' ! diffracted proton 

C— BEAM MOMENTA 

PBEAM1=27 . 6 

PBEAM2=820 . 

CALL HWIGIN 

M0DPDF(2)=-1 

iproc=19000 ! DIS hard subprocess 

NSTRU=6 ! selects HI pomeron structure routine 

Q2WWMN=l.E-6 ! minimum t 

Q2WWMX=4. ! maximum t 

YWWMIN=. 00299 ! minimum xp 

YWWMAX=. 00301 ! maximum x P 

Q2MIN=11.9 ! minimum Q 2 of DIS photon 

Q2MAX=12.1 ! maximum Q 2 of DIS photon 

The event record is shown in table 3. The diffracted proton is represented by entry 
5, labelled E-. Entry 4, labelled GAMMA is the pomeron. The incoming positron, 
entry 6, scatters off a UBAR quark in the pomeron, entry 7. The complete .F 2 D( ^ 
as generated by Pomwig with the default parameters is compared to the HI 
measurement Q in figure 2. Pomwig fails only in the highest (5 bins, that is at 
low diffractive masses, because the contribution from diffractively produced mesons 
is not included in Herwig . 

In PP collisions at the Tevatron, we consider two example processes, diffractive 
dijet production and Higgs production via pomeron-pomeron fusion. The single 
diffractive process PP — > PX is shown in table 4. For this process, HERWIG 
parameter iproc should be set to 11500. If diffractive W production were required, 
with the W decaying to electrons only, iproc would be set to 11451, and so on. 
The pomeron is represented in the event record by entry 4, labelled GAMMA, xp 
for this event is 0.0067. Entry 5 in the event record, labelled E+, is the scattered 
proton. In the hard subprocess, entry 6 is the parton from the pomeron, and entry 
7 is the parton from the proton. (3 for this event is 0.22. 

The double diffractive process PP — > PPH, for the Higgs decaying to bb is shown 
in table 5. iproc= 11605 selects Higgs production with the Higgs decaying only to 
bb. In this case, the two pomerons are entries 4 and 6 in the event record, labelled 
GAMMA. Entries 5 and 7 are the intact scattered protons, and entries 8 and 9 are 
the partons from the pomerons which enter the Higgs production hard subprocess. 
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Figure 2: The diffractive structure function F 2 as generated by Pomwig com- 
pared to the data of the HI Collaboration [[J. The blue dashed line shows the 
pomeron contribution, and the red solid line shows the sum of the pomeron and 
reggeon contributions. 
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EVENT 1: 27.60 GEV/C E+ ON 820.00 GEV/C E- PROCESS: 19000 
SEEDS: 313759 637835 STATUS: 100 ERROR: WEIGHT: 1.4238E-04 



—INITIAL STATE— 



IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


1 


E+ 


-11 


101 














0.00 


0.00 


27.6 


2 


E- 


11 


102 








4 


5 


0.00 


0.00 


-820.0 


3 


CMF 





103 


1 


4 








0.20 


0.41 


25.1 


4 


GAMMA 


22 


3 


2 











0.20 


0.41 


-2.5 


5 


E- 


11 


1 


2 











-0.20 


-0.41 


-817.5 



—HARD SUBPROCESS— 



IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


6 


E+ 


-11 


121 


8 


6 


11 


6 


0.00 


0.00 


27.6 


7 


UBAR 


-2 


122 


8 


10 


12 


10 


0.05 


0.10 


-0.6 


8 


HARD 





120 


6 


7 


9 


10 


0.07 


0.11 


27.0 


9 


E+ 


-11 


123 


8 


9 


14 


9 


3.06 


0.76 


22.8 


10 


UBAR 


-2 


124 


8 


7 


15 


7 


-3.01 


-0.66 


4.2 


11 


Z0/GAMA* 


23 


3 


6 


8 








-3.06 


-0.76 


4.8 



Table 3: Extract from a typical event record for diffractive DIS at HERA 



EVENT 1: 900.00 GEV/C E+ ON 900.00 GEV/C P PROCESS: 11500 
SEEDS: 313759 637835 STATUS: 100 ERROR: WEIGHT: 0.1482E+07 



—INITIAL STATE— 



IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


1 




E+ 


-11 


101 








4 


5 


0.00 


0.00 


900.00 


2 




P 


2212 


102 














0.00 


0.00 


-900.00 


3 




CMF 





103 


4 


2 








-0.18 


0.22 


-893.98 


4 




GAMMA 


22 


3 


1 











-0.18 


0.22 


6.02 


5 


E+ 


-11 


1 


1 











0.18 


-0.22 


893.98 


— Rl 


VRB SUBPROCESS— 


IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


6 




GLUON 


21 


121 


8 


7 


11 


9 


-0.04 


0.05 


1.35 


7 




GLUON 


21 


122 


8 


10 


15 


6 


0.00 


0.00 


-55.61 


8 




HARD 





120 


6 


7 


9 


10 


0.09 


0.03 


-54.25 


9 




GLUON 


21 


123 


8 


6 


20 


10 


4.12 


-0.99 


-2.04 


10 




GLUON 


21 


124 


8 


9 


26 


7 


-4.16 


1.04 


-52.22 



Table 4: Extract from a typical event record for the single diffractive process at 
Tevatron Run I 
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EVENT 1: 1000.00 GEV/C E+ ON 1000.00 GEV/C E- PROCESS: 11605 
SEEDS: 9876759 63576835 STATUS: 100 ERROR: WEIGHT: 0.5181E-08 



—INITIAL STATE— 



IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


1 


E+ 


-11 


101 








4 


5 


0.00 


0.00 


1000.00 


2 


E- 


11 


102 








6 


7 


0.00 


0.00 


-1000.00 


3 


CMF 





103 


4 


6 








0.42 


-0.24 


56.61 


4 


GAMMA 


22 


3 


1 











0.04 


0.17 


98.46 


5 


E+ 


-11 


1 


1 











-0.04 


-0.17 


901.54 


6 


GAMMA 


22 


3 


2 











0.38 


-0.41 


-41.85 


7 


E- 


11 


1 


2 











-0.38 


0.41 


-958.15 



—HARD SUBPROCESS— 



IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


8 
9 

10 


UBAR 
UQRK 
HIGGS 


-2 
2 

25 


121 
122 
120 


10 
10 

8 


9 
8 
9 


11 

13 

25 


9 

8 

25 


0.03 
0.33 
2.35 


0.15 

-0.36 

0.58 


82.80 

-36.53 

46.29 


— STRO 


NG HA 


DRON ] 


3ECA^ 


r S— 


IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


25 


HIGGS 


25 


195 


10 


25 


31 


32 


2.35 


0.58 


46.29 


— H/W/Z BOSON DECAYS— 


IHEP 


ID 


IDPDG 


1ST 


MOl 


M02 


DAI 


DA2 


P-X 


P-Y 


P-Z 


31 
32 


BQRK 
BBAR 


5 
-5 


123 
124 


25 
25 


32 
31 


33 
35 


32 
31 


-47.94 
50.29 


-1.24 
1.82 


48.97 
-2.69 



Table 5: Extract from a typical event record for Higgs production via pomeron- 
pomeron fusion at Tevatron Run II 



8 



References 

[1] G. Marchesini et al, Comp. Phys. Comm. 67 (1992) 465 ; G. Corcella et al, 
[hep-ph/9912396 

[2] http://www.pomwig.com 

[3] HI Collaboration: C.Adloff et al. Z.Phys. C74 (1997) 221. 
[4] J. F. Owens, Phys. Rev. D30 (1984) 943. 



9 



